Probing neutron correlations through nuclear break-up 
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The effect of initial correlations between nucleons on the nuclear break-up mechanism is studied. A 
quantum transport theory which extends standard mean-field approach is developed to incorporate 
short range pairing correlation as well as direct nucleon-nucleon collisions. A time evolution of the 
nuclear break-up from a correlated system leading to the emission of two particles to the continuum 
is performed. We show that initial correlations have strong influence on relative angles between 
particles emitted in coincidence. The present qualitative study indicates that nuclear break-up 
might be a tool to infer the residual interaction between nucleons in the nuclear medium. 
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Nuclei are self-bound systems formed of fermions in- 
teracting through the strong nuclear interaction. While 
many facets of nuclei could be understood in term of in- 
dependent particle motion, some aspects reveal internal 
correlations We consider here the so-called break-up 
process leading to the emission of nucleons to the contin- 
uum. Numerous dedicated models have been developed 
to account for this mechanism '^l. Among them, time 
dependent models based on the independent particle hy- 
pothesis have been shown to provide a good description of 
the nuclear as well as Coulomb break-up 01 ■ These ap- 
proaches, by neglecting two-body correlations could how- 
ever not provide appropriate theories when two nucleons 
are emitted from the same nucleus [1, H, 0] ■ Interferom- 
etry measurements are being now analyzed using rather 
schematic models and more elaborated theories, are 
highly desirable. 

The aim of the present work is twofold: (i) develop 
a microscopic quantum transport theory which incorpo- 
rates effects beyond mean-field like pairing correlations 
and/or direct nucleon-nucleon scattering in the medium, 
(ii) present a qualitative study of nuclear break-up and 
show that this mechanism can be a tool of choice for the 
study of correlations in nuclei. Similar challenges to (i) 
are being now addressed in strongly correlated electronic 
systems using Time-Dependent Density Functional The- 
ory (TDDFT) d, m. The Energy Density Functional 
(EDF) [l^ shares many aspects with DFT and is ex- 
pected to provide a universal treatment of static and dy- 
namical properties of nuclei [l0|, [U • Current EDFs start 
from an effective interaction (of Skyrme or Gogny type) 
to provide an energy functional, denoted £{p), where p is 
the one-body density matrix. Then, guided by the Hamil- 
tonian case, equations of motion are written in terms of 
the one-body density evolution given by ihdtp = [h[p],p], 
where h[p] = d£{p)/dp denotes the mean-field Hamil- 
tonian. To account for the richness of phenomena in 
nuclear dynamics (l2| . different extensions of mean- field 
have been proposed starting from the evolution: 



where denotes the effective vertex in the correlation 
channel, Tr2(.) is the partial trace on the second parti- 
cle. Ci2 denotes the two-body correlation defined from 
the two-body density pi2 as C12 = P12 — PiP2(l — ^12)- 
The indices refer to the particle on which the operator is 
applied (see for instance [HI, [l3|) while A12 is the per- 
mutation operator. Eq. ([T]) is generally complemented 
by the correlation evolution: 



ih 



dCvi 
dt 



[hi [p] + /i2 [p] , C12] + S12 + P12 + i?12 , (2) 
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[h[plp]+TT2[vl2,Cu 
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where again the indices in h refer to the particle to which 
the Hamiltonian is applied. Expression of B12 and P12 
and H12 , which can be found in [T^, [3] are guided by the 
BBGKY hierarchy [li]. These terms describe in-medium 
collisions, pairing and higher order effects respectively. 
When three-body correlations are neglected, the above 
theory reduces to the so-called Time Dependent Den- 
sity Matrix (TDDM) theory Coupled equations ^ 
and (21) have been directly applied to giant resonances in 
ref. iBl and more recently to fusion in ref. [l^. How- 
ever, applications are strongly constrained by the size 
of the two-body correlation matrix involved, fn addi- 
tion, similarly to ref. [17|, we encountered difficulties 
to obtain numerical convergence towards a stable corre- 
lated system. Therefore, an appropriate approximation 
should be made to render the TDDM theory more ver- 
satile. Keeping only B12 and projecting out the effect 
of correlation onto the one-body evolution leads to the 
so-called Extended TDHF theory with a non-Markovian 
collision term [H, Keeping only Pi 2 and assuming 
separable correlations leads to the Bogoliubov extension 
of the TDHF, i.e. TDHFB Both theories have been 
recently applied but require an [H, H^l • 

A different approximation is used here. We are inter- 
ested in nuclei at low excitation where pairing plays an 
important role. Similarly to TDHFB, we group single- 
particles into pairs, denoted by {a, a}, where \a) is ini- 
tially the time-reversed state of ja). We then assume that 
only components of and C12 between such pairs are 
different from zero. This approximation, called hereafter 
TDDM^, leads to important simplifications: (i) the num- 
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ber of correlation matrix components to be calculated is 
significantly reduced; (ii) the term H12 cancels out. In 
the basis where p is diagonal with occupation numbers 
Ua, i.e. p = |Q!)na(Q;|, the evolution reduces to 



ihdt\a) = h[p]\a) ; rio 



(3) 



where = (aa|wj2(l - ^i2)l/3^), C„/3 = (aa|Ci2|/?/3) 
and where the degeneracy of time-reversed states, i.e. 
Ua = TLa has been used. 

The TDDM^ incorporates correlations in the dynamics 
but also could be used to initialize a correlated system. 
A correlated nucleus is obtained in two steps. First the 
ev8 code ^21] is used to obtain single particle states which 
minimizes the Skyrmc EDF using SIII interaction. Sec- 
ond, the TDHF3D code 22], has been updated to incor- 
porate both the evolution of correlation (Eq. ^) as well 
as its coupling to the one-body density (Eq. U)). Equa- 
tions are integrated using a second order Runge-Kutta 
method where quantities on the right side of Eqs. (I3]l4|) 
are estimated at time t + At/2 to perform the evolution 
from t to t+lS.t. Note that this method insures the proper 
reorganization of the self-consistent mean-field when cor- 
relations built up. In particular, the mean-field polariza- 
tion due to correlations is accounted for. In this second 
step, following refs [H, [l3|, we make use of the Gell- 
Mann-Low adiabatic theorem [23^ and switch on adiabat- 
ically the residual interaction by v '^-^ jt ) = (l — e"*/'^) . 
The residual interaction is set to [10, : 



^12(^1:^2) = Wo (^1 - a |^p(i?)/poJ j S{fi - f2) (5) 

where R = {ri + r2)/2. The different parameters sets 
used here are given in table [J with po = 0.16 fm^'^. 



force 


Vo 


a 


/3 


Attractive 


-300 


1/2 


1 


Repulsive 


4-300 


1/2 


1 


Volume 


-159.6 







Surface 


-483.2 


1 


1 


Mixed 


-248.5 


1/2 


1 



TABLE I: Parameters of the different residual interactions 
used in this work. 

Applications are performed in a three dimensional 
Cartesian mesh of size (80 fni)'^ with a step of 0.8 fm 
and a time step of 0.45 fm/c. For r = 300 fm/c, a very 
good convergence of the adiabatic method, much better 
than for the full TDDM case has been achieved. 
For this first application of TDDM^, we consider the 
isotopic oxygen chain with an a core while correlations 



build up between neutrons belonging to the spd shells. 
After convergence, correlated systems have occupation 
numbers different from and 1 and non-zero correlation 
energy given by Econ — ^Tr(wj2C'i2)- Occupation num- 
bers are displayed in Fig. [TJ TDDM^ goes beyond pure 
pairing theory like HFB due to the inclusion of two par- 
ticles - two holes (2p-2h) terms. This is clearly illus- 
trated by the doubly magic ^^O nucleus where correla- 
tions do not cancel out and are in relatively good agree- 
ment with experimental observations of Refs. [1^. In 
table ini an estimate of the pairing gap defined as [l^] 
A = 2i?cori / J2a ~ '^a) Compared to the full 

TDDM and to the HFB cases. For the sake of compar- 
ison, the same core (^^O) and the same force as in ref. 
[13] has been used. We see that our results are glob- 
ally in agreement with the HFB and full TDDM results 
validating the approximation made in TDDM^ . 
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FIG. 1: Single-particle occupation numbers in oxygen iso- 
topes as a function of single particle energies. Lines represent 
Fermi functions fits. Open squares correspond to experimen- 
tal occupation probabilities in ^'^O [26i] . 





TDDM 113] 


TDDM^ 


HFB 


22q 


-3.1 MeV 


-3.5 MeV 


-3.3 MeV 


24 Q 


-2.7 MeV 


-3.1 MeV 


-3.4 MeV 



TABLE II: Effective pairing gaps for 22-24 Q deduced respec- 
tively from TDDM ^ and TDDM^ are compared to HFB 
results [11]. 

Correlated systems initialized with TDDM^ are then 
used to study the nuclear break-up leading to the emis- 
sion of two neutrons in coincidence. Intuitively, the two 
following scenarii have been proposed [2^: (i) if the two 
neutrons are initially close in position, both will feel the 
strong short range nuclear attraction of the reaction part- 
ner and will be emitted simultaneously at small relative 
angles (ii) if the two nucleons are far away in r-space, 
only one will undergo nuclear break-up. Then, the other 
nucleon might eventually be emitted isotropically from 
the daughter nucleus. Accordingly, large relative angles 
are expected between the two nucleons transmitted to 
the continuum in this sequential emission. To confirm 
this intuitive picture, the nuclear break-up dynamics is 
studied for an oxygen impinging on a ^"^Pb target. 
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FIG. 2; (Color online) One body density for three diflFerent 
steps of the dynamical evolution for an ^'^O +^°*Pb calcula- 
tion at 40 A.MeV. The circle represents the ^°*Pb projectile. 

The correlated nucleus is first initialized at the center 
of a 3D mesh of size (80 fni)'^. Concerning the dynamical 
step, the collision is simulated treating the collision part- 
ner as a one-body time-dependent external perturbation, 
by replacing first equation in ^ by: 

ihdt\a{t)) - {h[p{t)] + Vp{r,t)} \a{t)) (6) 

where Vp{f,t) stands for the projectile perturbation. 
Since we are considering here neutron emission. Coulomb 
effect is neglected and we assume that Vp is a mov- 
ing Woods-Saxon potential given by Vp{r,t) = Vo/(l + 
exp{(|r — ro{t)\ — i?pb)/a}) where ro(i) corresponds to 
the lead center of mass position. A simple straight line 
trajectory is used for r^it) corresponding to an impact 
parameter of 11 fm (grazing condition). The parameter 
i?Pb is set to 7.11 fm which corresponds to the ^"^Pb 
equivalent sharp radius while Vb=-50 MeV and a = 0.6 
fm. In the present application, we are interested in a 
specific reaction channel. If the full self-consistent calcu- 
lation is made, many other elastic and inelastic channels 
would be populated and would render the extraction of 
the relevant information difficult. To focus on the ef- 
fect of correlation on the break-up channel we assume 
that only correlated neutrons are affected by the colli- 
sions partner. Therefore, particles in the core as well 
as protons are fixed during the evolution. In addition, 
we assume that occupation numbers and correlation ma- 
trix elements are kept equal to their initial values. The 
calculation is performed in the oxygen frame and an ap- 
propriate transformation is made to simulate an oxygen 
beam and obtain the proper angular distributions in the 
laboratory frame. An example of one-body density evo- 
lution for the attractive force is shown in Fig. [2l In Fig. 
[21 three components can be distinguished: most of the 
neutrons remain in ^^O, part of the nucleons (inside the 
circles) have been transferred while the rest is emitted to 
continuum. The last component corresponds to nuclear 
break-up emission and has already been understood in [4] 
and observed experimentally [2^ . At the end of the evo- 
lution, nucleons from the inert core remaining in ^^O or 
nucleons transferred to the collision partner are removed 
from the following analysis. The relative angles be- 
tween two nucleons emitted in coincidence are then recon- 
structed from the correlation written in momentum space 
as: Ci2(pi,P2) = Ea/3'^"(Pi)'^a(P2)C„a^^(^/3(pi)0^(p2), 
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FIG. 3: Left: Relative angle correlation between neutrons 
at initial (top) and final (bottom) time of the evolution for 
an ^^O initialized with an attractive (full line) or repulsive 
(dashed line) residual interaction. Right: Initial (top) and 
final (bottom) relative angle correlation using the three dif- 
ferent residual interactions : a "Volume" (dashed line), "Sur- 
face" (full line) or "Mixed" (dashed-dotted line) residual in- 
teraction. All calculations are performed for an impact pa- 
rameter of fe = 11 fm. 



where the (jja denotes the fraction of the final time single- 
particle states emitted to the continuum. The distribu- 
tion of relative angles, C{9i2), is then deduced by sum- 
ming up contributions of all possible couples {pi,p2). 

To test the two scenarii, i.e. strong initial spatial cor- 
relation or anti-correlation, the "Attractive" and "Repul- 
sive" interaction have been used. The "Repulsive" inter- 
action is used to mimic the "cigar" like configurations, 
where the two correlated nucleons are well separated and 
on opposite side with respect to the core. Such a con- 
figuration is being now investigated experimentally using 
coulomb and/or nuclear break-up in lighter nuclei 6, 28]. 

The initial values of C{9i2) are given in top- left panel 
of Fig. [31 At initial time, the attractive force (close neu- 
trons in r-space) leads to large relative momentum (and 
therefore large relative angles). On the contrary, the re- 
pulsive force corresponds to small initial relative angles. 
In bottom- left panel, C'(6'i2) is displayed at the final time 
of the evolution. Several important remarks could be 
drawn. First (7(6*12) is largely modified compared to the 
initial ones underlying the importance of dynamical ef- 
fects. Therefore, experiments where two neutrons are 
detected in coincidence [Q| could only be used if proper 
transport model are developed. Second, we see that our 
calculation indeed confirms the intuitive picture. Strong 
initial correlation in space leads to small relative angle 
emission (solid line) while for initially well separated nu- 
cleons, relative angles are much larger. Comparison of 
our model results with experiments requires to sum up 
the contribution of different impact parameters. It is 
worth mentioning that the shape of the angular correla- 
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tion presented in Fig. [3] does not change much as the 
impact parameter or Wood-Saxon parameters are mod- 
ified. Therefore calculated cross-sections behave quali- 
tatively as in Fig. [3) In addition, experiments gener- 
ally present the ratio = P{Q\2) I P{P\)P{Q-2) where 
P(dv2) is the probability to emit two particles with a rel- 
ative angle C{Q\2) calculated in this article is the 
difference between the correlated emission and the inde- 
pendent emission. It can be compared directly to the 
experiment by normalizing it to the independent emis- 
sion and adding 1. In conclusion, experiments dedicated 
to the study of two nucleon emission to the continuum 
due to the nuclear break-up used in parallel with ded- 
icated transport models can provide a valuable method 
for the study of internal correlations in nuclei. 

We further investigated the sensitivity of relative an- 
gles with the initial correlations by comparing three real- 
istic residual interactions called "Volume" , "Surface" and 
"Mixed" in table U leading to the same scattering length 
(2^ . Corresponding initial and final angular correla- 



tions are respectively shown in top-right and bottom- 
right panel of Fig. [31 Although all forces correspond to 
the scenario (i), sizeable differences are observed in the 
amplitude of C(d\-2). Therefore, coincidence measure- 
ment of nucleon emission should be seriously considered 
in the near future as a tool to further constraint residual 
interaction used nowadays in EDF theories. It should 
however be kept in mind that a quantitative compari- 
son with experiments requires a careful analysis of the 
present model parameters which is underway. 
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